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Abstract. New data from HERA experiment on (diffractive) deep inelastic scattering has been used to 
parameterize nucleon and Pomeron structure functions. Within the Gribov theory, the parameterizations 
were employed to calculate gluon shadowing for various heavy ions and compared our results with predic- 
tions from other models. Calculations for d+Au collisions at forward rapidities at ultra-relativistic energies 
have been made and are compared to RHIC data on the nuclear modification factor. Results for gluon 
shadowing are also confronted with recent data on the nuclear modification factor at sfs = 17.3 GeV at 
various values of the Feynman variable xf, and the energy dependence of the effect is discussed. 

PACS. 12.40.Nn Regge theory, duality, absorptive/optical models - 13.60.Hb Total and inclusive cross 
sections (including deep inelastic processes) - 13.85.-t Hadron-induced high- and super-high-energy inter- 
actions (energy > 10 GeV) - 25.75.-q Relativistic heavy-ion collisions 



1 Introduction 

The fact that the nuclear structure function, F^ 4 , per num- 
ber of constituent nucleons is smaller than the structure 
function of a single nucleon, F% , for x < 0.1 1 is one of 
the most intriguing effects in modern high-energy nuclear 
physics. This effect is called nuclear shadowing. We define 
the nuclear ratio 



F 2 A 



(x,Q 2 ) 



(x,Q 2 



(1) 



which is smaller than unity in the shadowing region. The 
nuclear ratio has been measured for many nuclei and re- 
veals also an interesting structure for x > 0.1 (see pQ for 
further references). In high-energy hadron-nucleus colli- 
sions we probe the partonic structure of the nucleus at 
low values of x. For a given energy the smallest available 
values of x of the partons in the nucleus are to be found 
in the fragmentation region, where xp = p z /p™ ax of the 
observed particle is close to 1, which corresponds to large 
pseudorapidity 77. 

An understanding of the so-called cold nuclear effects, 
or initial state effects, in hadron-nucleus collisions are there- 
fore an important benchmark for nucleus-nucleus colli- 
sions. 
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1 In the infinite momentum frame, x is the longitudinal mo- 
mentum fraction of a parton in the nucleon 



A significant change in the underlying dynamics of a 
hadron-nucleus collision takes place with growing energy 
of the incoming particles. At low energies, the total cross 
section is well described within the probabilistic Glauber 
model [5], which only takes into account elastic rescat- 
tcrings of the initial hadron on the various nucleons of 
the nucleus. Elastic scattering is described by Pomeron 
exchange. At higher energies, E > E cr n ^ tyin [iRa corre- 
sponding to a coherence length 



1 



2 TOjV x ' 



(2) 



the typical hadronic fluctuation length can become of the 
order of, or even bigger than, the nuclear radius, Ra, and 
there will be coherent interaction of constituents of the 
hadron with several nucleons of the nucleus. The sum of all 
diagrams was calculated by Gribov . In this framework, 
the diffractive intermediate states has to be accounted for 
in the sum over subsequent rescatterings. The space-time 
picture analogy to the Glauber series is lost, as the inter- 
actions with different nucleons of the nucleus happens in- 
stantaneously. The phenomenon of coherent multiple scat- 
tering is referred to as shadowing corrections. 

An additional effect which comes into play at high 
energies, is the possibility of interactions between soft 
partons of the different nucleons in the nucleus. In the 
Glauber-Gribov model this corresponds to multi-Pomeron 
interactions. These diagrams are called enhanced diagrams 
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0j, and can also be understood as interactions between 
strings formed in the collision. E.g. the triple-Pomeron 
vertex is proportional to A 1 / 3, in hA collisions, and so it 
becomes very important for collisions on very heavy nu- 
clei. 

In what follows we will describe a model for hadron- 
nucleus collisions in Section[5] In Section|3we will discuss 
the experimental data on both inclusive and diffractive 
deep inelastic scattering (DIS) cross-sections measured at 
HERA used as input to the model. In Section 01 we will 
present the results for the shadowing ratio, and also com- 
pare them to recent experiments at the Relativistic Heavy 
Ion Collider (RHIC) in Sectional In SectionElwe compare 
our calculations to recent measurements done at much 
lower energies, y/s — 17.3 GeV, and discuss the energy 
dependence of the shadowing effect. We summarize and 
conclude in Section 



convenient as it guarantees the disappearance of nuclear 
shadowing at x ~ 0.1 as in experimental data. Coherence 
effects are taken into account in the form factor 



F A (t r , 



d bJ Q (b\/~t mm ) T A (b) 



(5) 



which is equal to 1 at x — > and decreases with increasing 
x due to the loss of coherence for x > (2m ^R a) _1 , see 

(here we have put t m i n = —m 2 N xp). 

The second order elastic cross section in (@J is obvi- 
ously negative and will lead to a reduction of the total 
cross section. In the small x region, it is also necessary to 
include higher orders terms in © in order not to violate 
unitarity of the total cross section, as was noted in 

Summation of all terms in @ is model dependent. The 
Schwimmer unitarization 1D| for the total hA cross sec- 
tion, which also sums up all Pomeron tree diagrams, is 
used to obtain 



2 The Model 

We will assume that the nucleus consists of A independent 
nucleons, in the spirit of the Glauber model. The scatter- 
ing amplitude of an incoming hadron on a nuclear target, 
can be expanded in a multiple scattering series |D 



a a = Aa 



N 



r (2) 



(3) 



where the first term is simply the Glauber elastic rescat- 
tering (see Fig. The second term in © is related to 
diffractive DIS through the AGK cutting rules (Sj. It is 
given by |Zj 



r (2) _ 



-4ttA(A-1) J d 2 bTl(b) 



AM 2 



M- 



dM 2 dt 



F A (t„ 



(4) 



where M 2 is the mass of the diffractively produced inter- 
mediate state, is the normalized nuclear thickness 
function and daJ^ N /dM 2 dt is the differential cross section 
for diffractive dissociation of the hadron. In the second 
integral, M^ ax — Q 2 (xp ax /x — 1) is found by demanding 
a large rapidity gap in the diffractive dissociation. Cal- 



culations are made both for 
for x% ax 



bp 



= 0.1 as in [5] and 



T Sch 

7*A 



A< 



7* N 



d 2 b 



T A (b) 



1 + [A-l)f{x,Q 2 )T A {b) 



(6) 



where f(x,Q 2 ) is the effective shadowing function. Fol- 
lowing jl 1U12| in choice of parameters and factorization, 
one can get the shadowing function as 



f(x,Q 2 ) = 4tt ff a *dx P B(xip) 



F^( X1P ,Q 2 ,(3) 



X F ^ (trnir! 



F 2 (x,Q 2 
)■ 



(7) 



Here B( X]P ) = 0.184- 0.02 In (xjp) fm 2 . The lhs. in © is 
defined as the shadowing ratio R h (A/N) (x). 

The structure function F 2 and the diffractive structure 
function i^p of the single nucleon are taken as input from 
experiment. The extension to the nuclear case is therefore 
parameter-free except for the unitarity constraints leading 
to eq. JBJ). This is a remarkable feature of the Glauber- 
Gribov model 1 1 H j . 

The model is valid for low values of x < 0.01 and 
intermediate Q 2 < 10 GeV 2 . In what follows, we will ne- 
glect the Q 2 dependence of the model, knowing that the 
presence of a strong Q 2 -dependent term is not required to 
describe nuclear data at low Q 2 which is relevant for our 
present considerations 



0.03 as in [Hj, although the former is more 3 | nc | us j ve and diffractive data 




Fig. 1. The single and double scattering from to the total 
cross section (the diffractive intermediate states are on mass 
shell). 



In deep inelastic scattering (DIS) the structure function 
of a nucleon is related to the total cross section of 7* TV in- 
teraction through factorization at high scales valid in per- 
turbative QCD. The structure function holds information 
about the partonic content of the nucleon, and is given by 
a sum of parton distribution functions (PDFs) 



F2 (x, Q 2 ) 



oc 



E 

=g,u,d,s. 



(8) 



where the sum is over all types of partons. Similar to 
the inclusive DIS case, a factorization theorem has been 
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proved in perturbative QCD to hold for diffractive struc- 
ture functions fTl] . The diffractive structure function 
in eq. J7J) is given by 



F% ( X]P ,Q 2 ,I3) = fjpixjp) Ff (f3,Q 2 



(9) 



where fp is the i-integrated Pomeron flux and we have 
assumed so-called Regge factorization. The Pomeron flux 
factor is defined as 



At 



2a F (t)-l 



dt 



(10) 



where we assume a linear Pomeron trajectory, ap(t) = 
ap(0) + a'pt. The values of the parameters are taken to 
be ajp(0) = 1.173 and a'p = 0.26 GeV~ 2 , and we put 
B = 4.6 GeV~ 2 (see |H| for details). Ff is called the 
structure function of the Pomeron, and so it is possible to 
introduce a partonic structure of the Pomeron |16j . as in 
eq. ®. 

For both inclusive and diffractive DIS, the gluon con- 
tent clearly dominates over the quark one at low x and 
intermediate Q 2 values relevant to our present considera- 
tions. In what follows we will only consider the gluon par- 
ton distribution functions of the nucleon and Pomeron. 
Quark shadowing was discussed in 0. The ratio under 
the integral in eq. can, in other words be understood 
as the density of gluons in the Pomeron compared to the 
density of gluons in the nucleon! 

For the nucleon structure function we will use the next 
to leading order (NLO) ZEUS-S QCD fit of the gluon PDF 
[T7] at Q 2 — 7 GeV 2 . The gluon diffractive parton distri- 
bution function (dPDF) was measured for intermediate 
Q 2 at both HERA experiments: ZEUS and HI. The sit- 
uation regarding the gluon content of the Pomeron is at 
present rather uncertain with a large discrepancy between 
the presented results. The origin of the discrepancy is un- 
known. We will make calculations for three different fits 
which we will briefly mention 

- HI NLO QCD 2002 fit Q3] (large-rapidity-gap data, 
gives also a good description of heavy quark and dijet 
production) 



- H1 NLO QCD 2002 (LRG) 




" H1 & ZEUS 2005 (Newman 


Gt al.) 




" H1 NLO QCD 2006 (FPS) ..**' 
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Fig. 2. Nuclear shadowing ratio i?(Pb/N) for different param- 
eterizations of the diffractive gluon parton distribution func- 
tion. 



- HI & ZEUS 2005 fit [IB] (HI NLO QCD analysis of 
ZEUS Mx-data) 

- HI NLO QCD 2006 fit [T31I27)] (forward proton spec- 
trometer data, fit A was used) 

The gluon dPDF for the Pomeron was parameterized at 
fixed Q 2 = 6.5 GeV 2 (Q 2 = 8.5 GeV 2 in the second case). 
Both inclusive and diffractive distributions were fitted by 
a simple function 

xf g {x,Q 2 ) = xf g (x) = Ax'^l-x) 1 , (11) 

where A, 5 and 7 are fitting parameters. 



4 Nuclear shadowing ratio 

Numerical calculations of the nuclear shadowing ratio for 
Pb, defined in eq. ©, for different parameterizations of 
the gluon dPDF as described in the previous section are 
presented in Fig. FJ The fit to ZEUS data [H] predicts 
weakest gluon shadowing, while the fits to the both of the 
HI datasets are compatible with each other and predict 
a much stronger shadowing effect. The strongest gluon 
shadowing is obtained for the HI NLO QCD 2002 fit [E], 
which is almost twice as big as the ZEUS one flS] for the 
whole range of x. 

Gluon shadowing for various heavy ions (Ca, Pd and 
Pb) calculated with © is presented in Fig. [3] (HI NLO 
QCD 2002 fit [El is used). The effect is strong at small 
x, and disappearing at x = Xp ax . This is a consequence 
of the coherence length in the form factor (JSJ), and the 
vanishing integration domain in 0. Gluon shadowing is 
as low as 0.2 for the Pb/N ratio at x ~ 10~ 5 . 

A comparison of our results for Pb/nucleon ratio at 
Q 2 = 6.5 GeV 2 with x^ ax = 0.03, with other models, cal- 
culated at Q 2 = 5 GeV 2 , is presented in Fig. H (HI NLO 
QCD 2002 fit [12] is used). In [21], the authors have cal- 
culated shadowing within the BFKL formalism [22], while 
|23j is based on a parameterization of pp-data. The au- 
thors of [H] make their calculations within a similar frame- 
work as the presented model. 

For x < 10~ 3 our model predicts stronger gluon shad- 
owing compared to [21] (dashed-dotted line) and [H] (dot- 
ted line), while [23] (dashed line) predicts the strongest 
effect down to x ~ 10~ 4 . Our calculations are close to the 
predictions of [H] for x > 10 -3 for this choice of xfi 
while we are consequently below the predictions of |21| . 



y,max 
hp 



5 Shadowing effects in d+Au collisions 

The model is now employed to study particle production 
in d+Au collisions at RHIC energy s/s = 200 GeV. There 
has been observed an increasing suppression of the nuclear 
modification factor (NMF) 



RdAu — 



1 d z N dAu /dprdrj 



(Neon) tfNZi/dprdfl 



(12) 
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• Pb/N 
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Fig. 3. Gluon shadowing for heavy ions. Closed (open) sym- 
bols are for x£ ax =0.1 (0.03). 



2 1.4 
I 12 
1 

0.8 
0.6 
0.4 
0.2 



— Armesto et.al: BFKL ladders 

— New HIJING 

Frankfurt et. al: Leading-twist shadowing 

— Glauber-Gribov, jj"" = 0.03 



10 4 



Fig. 4. Comparison of our results the Pb/nucleon ratio with 
other models, at fixed Q 2 . 



with increasing pseudorapidity of the observed particle 
|2*4*| . Various models have been utilized to explain the ob- 
served features, the most successful of which is the Color- 
Glass- Condensate (CGC) model which assumes gluon 
saturation for the kinematical domain reached in hadron- 
nuclcus collisions at RHIC. It is instructive to point out, 
that the model presented here does not assume a satura- 
tion scale, like in [53]. Yet, as we have already mentioned, 
Pomeron interactions are taken into account to preserve 
unitarity of the scattering amplitude. In this sense it should 
give similar results as models assuming gluon ladder fu- 
sion. 

In the Glauber-Gribov model, the multiplicity reduc- 
tion due to shadowing compared to the simple Glauber 
model is easily obtained in a factorized form. The theo- 
retical prediction is given by |26) 



jjtheo 
n dAu 



Rr n (x p )R s A l h (x t ) , 



(13) 



where the deuteron will be treated as a point particle in 
impact parameter space, but with the shadowing found 
from ©• The collision is described by the following jet 
kinematics 



Cp(i) 



cpx e 



(14) 



for the projectile (target) cc-value respectively. In eq. i|14[l 
Pt is the transverse momentum of the particle, and we 
assume that most of the high-pr particles come from jets 
c times more energetic. 
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1.2- 
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0.6 
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c = 5 



i* * 
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R d Au(n=o)/R d h r u (n=o) 
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1 2 3 4 5 6 

p T , GeV/c 

Fig. 5. Shadowing effects for d+Au collisions \/s = 200 GeV 
at midrapidity for c = 5. 



An important and well-known effect that is not taken 
into account in the model presented here, is the Cronin 
effect [27112%] . or j>T-broadening of the produced particles, 
which leads to an effective enhancement of the NMF seen 
at midrapidity for px > 2 GeV 24 29 . In what follows we 
will assume that this effect stays the same for all pseudo- 
rapidities. 

We therefore extract the gluon shadowing effects in 



the NMF at r, = by defining Bffi? = [KaJ Ka°u\ v ^ 
shown in Fig.Elfor c = 5. As expected, the shadowing dies 
out at high px ■ 

The multiplicity reduction merely due to shadowing 
will then appear as we compare the NMF at forward ra- 
pidities, T] = 1, 2.2, 3.2 to R^ZT- This can be quantified 
by the double ratio 



R 



[RdAujri 

r>norm 
K-dAu 



(15) 



where the NMF in the numerator is taken at a constant 
rapidity slice. 

The double ratio, as defined in (|15fl . for the NMF taken 
from [23] is plotted together with the predictions of the 
Glauber-Gribov model from eq. l|13f) in figs.|H|and[7Jfor two 
different values of the parameter c. Calculations are made 
for the three gluon dPDF parameterizations described in 
Section Statistical errors are denoted by the thick solid 
line, while the systematic and statistical errors added up 
quadratically are denoted by the dashed line. The choice 
of c does not seem to affect the result. 

There is good agreement with experimental data for 
all the gluon dPDF parameterizations. The HI NLO QCD 
2002 fit seems to be most consistent with the data, and 
the agreement is better at higher values of pseudorapidity. 
Although the ZEUS fit (HI & ZEUS 2005) is almost a 
factor of two smaller for the whole range of possible x, the 
agreement is surprisingly good. This gives evidence for a 
dominating shadowing contribution to the suppression of 
the NMF at forward rapidities at RHIC energies. 
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Fig. 6. NMF ratio, c = 3 



Fig. 7. NMF ratio, c = 5 



6 Nuclear shadowing at SPS 



We have also calculated the nuclear shadowing ratio for 
lower energies, namely at maximal SPS energy yfs — 17.3 
GeV. In Fig. [H] we show the calculation of the Glauber- 
Gribov model compared to data on charged pion produc- 
tion for fixed xf — 0.375 taken from [3U| for the double 
ratio defined in eq. {H) (the HI NLO QCD 2002 fit was 
used) . The curves are for two values of the parameter c, 
and we see that the shadowing disappears quickly with 
increasing c. 

Obviously, the effect of gluon shadowing is not suffi- 
cient to explain the observed suppression at this energy. 
The solid curve in Fig. [H] can be taken as the maximum 
value of the effect. The suppression in the experimental 
data is strongest for small pr- 

Theoretical considerations on the reasons for the sup- 
pression in the NMF at SPS energies are out of the scope 
of the present paper, but will be followed up in the near- 
est future An extremely interesting fact is that the 
suppression at SPS is almost of the same magnitude as at 
RHIC energies. Since gluon shadowing is expected to be- 
come significant with growing energy of the reaction, there 
is apparently another mechanism present which is respon- 
sible for the suppression. This mechanism is related to 
the energy-scale relevant for coherent scattering; at these 
energies a large fraction of the Fock-state of the incom- 
ing hadron will rescatter as in the Glauber model |31) . 
Energy-momentum conservation effects, which violate the 
AGK cutting rules, will play a dominant role |33|. 
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Fig. 8. Shadowing ratio for y/s — 17.3GeV in the fragmenta- 
tion region. 

7 Conclusions 

We have presented results for the nuclear shadowing ra- 
tio, eq. ©, and particle production in hadron- nucleus col- 
lisions for various parameterizations of the gluon dPDF 
within the Glauber-Gribov model, and compared them 
to experimental data on the nuclear modification factor 
measured at two different energies, yfs = 17.3 GeV and 
yfs — 200 GeV. Our calculations of the nuclear shadow- 
ing ratio is consistent with other models predicting large 
gluon shadowing for low values of x and intermediate Q 2 . 

For RHIC energies, the suppression of the NMF ob- 
served at non-zero values of pseudorapidity is well de- 
scribed by gluon shadowing within the Glauber-Gribov 
model. The experimental data on gluon dPDF is still quite 
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uncertain, and introduces a large spread in the theoreti- 
cal prediction. The agreement with experimental data is 
reasonable for all the presented parameterizations. 

We would like to underline that the results we have 
presented can be viewed as an upper bound of the effect 
of gluon shadowing in hadron- nucleus collisions. The au- 
thors of |34) have done calculations within a similar frame- 
work with a different choice of kinematics than in eq. (|14(l , 
resulting in a much weaker shadowing effect. This discrep- 
ancy is important to resolve in the nearest future. 

In |34| there is also an important remark on the ex- 
perimental data from BRAHMS [23] at the two most for- 
ward rapidities. The fact that only negative particles, h~ , 
are measured leads effectively to an enhancement of the 
NMF because of isospin effects. In order to compare to the 
correct NMF, one should reduce the experimental data at 
r\ = 2.2 and 3.2 by a factor of ~ 2/3. This will lead to a less 
impressive agreement of our model, yet it will not change 
the conclusions regarding gluon shadowing at RHIC. 

At SPS, the gluon shadowing is not responsible for 
more than 10 % of the total suppression. The important 
fact, however, is the observed large suppression up to pt = 
2 GeV/c at such low energies. The suppression seems in 
fact to be approximately of the same magnitude as at 
RHIC, where y/s is a factor often larger. The suppression 
is caused by energy-momentum conservation which should 
stay constant with growing energy. This may indicate that 
our estimates of gluon shadowing are too large. 

The energy-dependence of the suppression is related to 
the underlying space-time dynamics of the collision, and is 
therefore a crucial test for theoretical models. With new 
low-energy data on NMF in the forward region a com- 
parison of the effect can now also be done for xf > 0. 
This gives an opportunity to study the interplay of differ- 
ent effects leading to suppression/enhancement of particle 
spectra in much more detail. 

The presented model can also be used to calculate the 
expected suppression in heavy-ion collisions at LHC ener- 
gies 122]. 
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